PERFORMANCE  RANKING  OF  FINNED 
TUBULAR  HEAT  EXCHANGER  SURFACES 


James  Ernest  Baskervi lie 


MRUMMf 
SMMUHE 


PERFORMANCE  RANKING  OF 
FINNED  TUBULAR  HEAT  EXCHANGER  SURFACES 

by 

JAMES  ERNEST  BASKERVILLE 

// 


B.S.,  UNITED  STATES  NAVAL  ACADEMY 
Annapolis,  Maryland 
(1969) 


SUBMITTED  IN  PARTIAL  FULFILLMENT  OF  THE  REQUIREMENTS 
FOR  THE  DEGREE  OF  OCEAN  ENGINEER  AND  MASTER  OF  SCIENCE 

IN  MECHANICAL  ENGINEERING 

at  the 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 

May  1977 


"■W  POSTGRADUATE  SCNOOt 


PERFORMANCE    RANKING   OF 

FINNED   TUBULAR   HEAT    EXCHANGER   SURFACES 

by 

JAMES  ERNEST  BASKERVILLE 

Submitted  to  the  Department  of  Ocean  Engineering  on  May  12, 
1977  in  partial  fulfillment  of  the  requirements  for  the  degree 
of  Ocean  Engineer  and  Master  of  Science  in  Mechanical  Engin- 
eering. 

ABSTRACT 

The  traditional  method  of  data  presentation  for  heat  ex- 
changer surfaces  does  not  permit  comparison  of  individual  sur- 
face types  in  any  simple  manner.   This  data  is  most  commonly 
presented  in  terms  of  heat  transfer  coefficients  and  friction 
factors  referenced  to  the  exposed  area  as  a  function  of  Rey- 
nolds number  based  on  the  minimum  free  flow  area. 

Soland  [5]  proposed  a  method  of  surface  comparison  for 
plate-finned  heat  exchanger  surfaces  in  which  the  heat  trans- 
fer coefficient  and  friction  factor  is  referenced  to  the  base 
area  and  Reynolds  number  is  based  on  the  open  flow  area,  as 
though  the  enhanced  surfaces  were  not  present. 

Theory  and  principle  inherent  in  the  Soland  proposal  are 
used  to  develop  a  method  of  surface  comparison  for  finned  tube 
surfaces  as  presented  in  (a)  Kays  and  London  [2] ,  and  (b)  sev- 
eral supplemental  surfaces  furnished  by  the  Trane  Corporation. 
The  derived  comparison  method  is  evaluated  by  application  to 
a  practical  crossflow  finned  tubular  heat  exchanger  design. 

Appendix  III  provides  a  method  for  sizing  crossflow 
finned  tubular  heat  exchangers. 
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NOMENCLATURE 


a.        Minimum  free  flow  dimension  between  tubes  measured 
d 

diagonally  between  adjacent  rows  from  base  to  base 
as  if  fins  were  not  present;  see  Figure  4,  [FT] 

a.         Minimum  free  flow  dimension  between  tubes  measured 
transversely  to  the  direction  of  flow  from  base  to 
base  as  if  fins  were  not  present;  see  Figure  4, 
[FT] 


*b 


Heat  transfer  area  of  the  bare  base  tubular  sur- 
faces; equals  length  times  heated  perimeter  times 

2 
number  of  tubes;  [FT  ] 


A         Exchanger  minimum  free  flow  area,  defined  by  (la) , 

c  2 

[FT  ] 


Af        Constant  from  linearized  equation  for  friction  fac- 

B 
tor  (f)  of  the  form:  A  RE 


A_        Exchanger  flow  passage  frontal  area  ignoring  any 

2 
enhancing  surfaces;  see  Figure  4,  [FT  ] 

2 
A, /A      Ratio  of  fin  area  to  total  transfer  area,  [FT  ] 

2 
Af        Exchanger  total  frontal  area,  [FT  ] 

A.         Constant  from  linearized  equation  for  Colburn  mod- 

B 

ulus  (j)  of  the  form:  A  RE 


A„        Total  heat  transfer  area  of  exchanger  on  one  side, 

2 
defined  by  3  times  total  volume,  [FT  ] 


b         Longitudinal  center-to-center  tube  spacing,  see 
Figure  4,  [FT] 

C         Flow  stream  capacity  rate  [WC  ] ,  [BTU/HR-°F] 

Cn        Constant  used  in  the  computation  of  friction  factor 
(f)  for  flow  over  smooth  tube  banks;  defined  by 
(32) 


C.        Constant  required  for  the  computation  of  radial 
fin  efficiency;  defined  in  Figure  5. 


C         Specific  heat  at  constant  pressure,  [BTU/LBM-°F] 

XT 

D.        Constant  required  for  the  computation  of  radial 
fin  efficiency;  defined  in  Figure  5. 


D.         Hydraulic  diameter,  defined  by  (la) ,  [FT] 

D„        Nominal  diameter,  defined  by  (lb) ,  [FT] 
N 

F         Transverse  center-to-center  tube  spacing;  see  Fig- 
ure 4,  [FT] 


friction  factor  based  on  total  area  (A  ) ;  defined 
by  (4a) 


fXT        friction  factor  based  on  base  area  [A,  ]  ;  defined 
N  D 

by  (4b) 


f         friction  factor  based  on  total  area  (A_)  for  a 

S  T 

smooth  surface,  defined  by  (31)  or  (32) 
g0         32.174  [LBM/LBF]  [FT/SEC2] 


G         Mass  flux  based  on  minimum  free  flow  area,  [A  ] , 
c  ^  c 

defined  by  (2a),  [LBM/HR-FT  ] 


G         Mass  flux  based  on  free  flow  area,  [A_] ,  defined 

N  —  r 

by  (2b),  [LBM/HR-FT^] 


Heat  transfer  coefficient  based  on  total  area,  [Ar 


T 
defined  by  (5a),  [BTU/HR-FT^- °F] 


.2  o 


h-.        Heat  transfer  coefficient  based  on  base  area,  [A,] ; 
defined  by  (5b),  [BTU/HR-FT  -°F] 


Colburn  j-f actor  based  on  total  area,  [A  ] ;  define 
by  (7a) 


JN        Colburn  j-f actor  based  on  base  area,  [A  ] ;  defined 
by  (7b) 

j         Colburn  j-f actor  based  on  total  area,  [A  ] ,  for  a 

S  x 

smooth  surface,  defined  by  (28) 

K         Thermal  conductivity  [ BTU/HR-FT- °F] 

I  Fin  height  defined  by  r  -r  ;  for  rectangular  fin 

sheets,  see  Figure  5,  [FT] 

L         Flow  length  along  the  axis  of  the  tubes  (i.e.  for 
the  fluid  inside  the  tubes) ;  for  fluid  flowing 
normal  to  the  tube  banks,  defined  as  N  times  b; 
[FT) 

m         Component  of  fin  efficiency  (n ,   L) ;  defined  by  [10] 

N         Number  of  tubes 


N_        Symbol  used  in  reference  [2]  for  RE 


Nu        Nusselt  number;  defined  by  (6a) 


Nu„       Nusselt  number;  defined  by  (6b) 


NTU       Number  of  transfer  units;  defined  by  (21) 


H/  L      Fi-n  temperature  effectiveness  of  a  longitudinal 
fin  of  rectangular  profile;  defined  by  (9) 


An„„_    Correction  factor  applied  to  t),    (9)  to  permit  the 
corr  q— jj 

use  of  the  simplier  [Tanh  ml] /ml  equation  for  fin 
temperature  effectiveness  in  lieu  of  Bessel  Func- 
tion type  equations. 


n         Total  surface  temperature  effectiveness,  defined 
by  (8) 


Pr        Prandtl  Number  [yC  /K] 

2 

P         Pressure  [LBF/FT  ] 


Pumping  Power  [HP] 


2 
APf       Friction  pressure  drop  [LBF/FT  ] 


q         Heat  transfer  rate  [BTU/HR] 


q/A        Heat  Flux  [BTU/HR- FT2] 


RE        Reynolds  number  based  on  minimum  free  flow  area, 
[Ac] ,  defined  by  (3a),  Labeled  NR  by  ref  [2] 
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REN       Reynolds  number  based  on  free  flow  area,  [A  ] ,  de- 
fined by  (3b) 

r         Radius  of  enhanced  (finned)  surface,  [FT] 

r         Outside  radius  of  base  tube,  [FT] 

r.        Hydraulic  radius,  equals  hydraulic  diameter  divided 
by  four  for  a  round  tube,  [FT] 

T  Temperature,  [°F] 

U  Overall  heat  transfer  coefficient  [BTU/HR-FT  -°F] 

V  Volume  [FT3] 

W  Mass  flow  rate,  [LBM/HR] 

X,Y,Z  Principle  dimensions  of  heat  exchanger,  [FT] 


X_        Ratio  of  diagonal  tube  spacing  to  tube  outside 
diameter 


X_        Ratio  of  longitudinal  tube  spacing  to  tube  outside 
diameter,  equals  b  divided  by  D 

X_        Ratio  of  transverse  tube  spacing  to  tube  outside 
diameter,  equals  F  divided  by  D 


Miscellaneous 

a        In  tube  geometry,  angle  defined  by  the  longitudinal 
centerline  and  the  diagonal  center line;  see  Figure 
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7,  [degrees] 


Ratio  of  total  transfer  area  on  one  side  of  the 

exchanger  to  total  volume  of  the  exchanger,  label- 

2    3 
ed  as  a  in  reference  [2] ,  [FT  /FT  ] 


A,  -0.00259,  Constant  used  to  calculate  CA 

1  A 

An  -0.00840,  Constant  used  to  calculate  C. 

2  A 

A-  -0.02426,  Constant  used  to  calculate  D. 

A.  -0.04521,  Constant  used  to  calculate  D. 

6  Fin  thickness,  or  average  fin  thickness  if  tapered, 
[FT] 

e  Exchanger  effectiveness 

U  Viscosity  [LBM/HR-FT] 

p  Density  [LBM/FT3] 


Subscripts 

d         diagonal 

e         enhanced  surface 

m         heat  exchanger  fin  material 

0         outside 
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Subscripts  (continued) 
os        one  side 
s         smooth  surface 
t         transverse 


Total 

Finned  side,  colder  fluid  flowing  over  finned  tube 
banks 

Tube  side;  hotter  fluid  flowing  inside  smooth  tubes 
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I.   INTRODUCTION 

A.   Purpose 

The  objective  of  any  heat  exchanger  comparison  method  is 
to  enable  the  designer  to  select,  from  an  often  overwhelming 
number  of  enhanced  surfaces,  that  unique  surface  which  can  be 
termed  "optimum".   The  comparison  technique  should  be  in  the 
form  of  a  logical,  easily  implemented,  inherently  accurate 
procedure. 

It  should  be  recognized  from  the  onset  that  regardless 
of  the  definition  of  "optimum"  selected,  the  surface  labeled 
as  such  may  not,  in  fact,  be  the  best  for  a  given  practical 
situation.   Space  limitations  or  cost  restraints  frequently 
restrict  the  designers  flexibility.   While  a  particular  surf- 
ace may  be  "optimum"  from  a  least  volume  criteria,  for  exam- 
ple, it  may  well  exceed  the  cost  limitations  or  require  an 
unacceptably  long  dimension.   Violation  of  either  of  these 
practical  constraints  could  well  force  the  selection  of  an 
alternate  surface. 

Soland  [5]  proposed  a  method  of  comparison  that  permits 
performance  comparisons  of  all  of  the  Kays  and  London  [2] 
plate  finned  type  heat  exchanger  surfaces  on  four  different 
basest 

a.  Same  exchanger  shape  and  volume 

b.  Same  exchanger  volume  and  pumping  power 

c.  Same  pumping  power  and  NTU 
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d.  Same  volume  and  NTU 

Sheldon  [6],  (a)  applied  Soland's  method  to  practical 
heat  exchanger  design  problems;  (b)  evalueated  the  method  ef- 
fectiveness; and  (c)  performed  a  comparative  analysis  on  plate 
finned  surfaces  in  addition  to  those  found  in  Kays  and  London. 

The  purpose  of  this  paper  is  to: 

1.  adapt  Soland's  method  for  application  to  crossflow 
finned  tube  type  surfaces  presented  in  (a)  Kays  and  London; 
and,  (b)  several  supplemental  surfaces  furnished  by  Trane 
Corporation; 

2.  evaluate  the  method  effectiveness  by  application  to 
practical  crossflow  finned  tubular  heat  exchanger  design 
problems. 

B.  Background 

Crossflow  finned  tubular  heat  exchangers  have  found  app- 
lication as  gas  turbine  plant  intercoolers,  aircraft  engine 
and  electronics  coolers,  air  conditioning  units,  heat  ex- 
changers, and  numerous  other  uses.   These  applications  most 
often  involve  gas-liquid  service.  [See  Figure  1] 

The  performance  of  a  given  heat  exchanger  in  which  the 
ratio  of  the  heat  transfer  coefficients  between  the  two  work- 
ing fluids  is  greater  than  about  three  is  markedly  improved 
by  the  employemnt  of  enhanced  surfaces.   The  increased  sur- 
face area  inherent  with  the  employment  of  fins  yields  a  net 
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FIGURE  1;   Finned  Tubular  Heat 
Exchangers 
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improvement  in  heat  transfer  despite  lower  heat  transfer  co- 
efficients.  Finned  tubular  heat  exchangers  have  surface-to- 

2    3 
volume  ratios  (0)  ranging  between  2  0-200  ft  /ft  . 

Reference  [2]  will  be  considered  the  primary  source  for 
data  related  to  the  design  and  testing  of  finned  tubular  sur- 
faces, while  supplemental  data  has  been  provided  by  the  Trane 
Corporation. 

Optimum  design  is  often  used  synonomously  with  "lowest 
cost".   Cost  is  frequently  defined  as  (a)  amortized  acquisi- 
sition  cost  related  to  the  initial  design  and  fabrication  of 
the  component,  (b)  operating  costs,  defined  as  amortized  acq- 
uisition cost  plus  operating  expenses  and  (c)  cost  to  the 
system  in  terms  of  either  weight  or  volume. 

Consider  aircraft  engine  component  design,  for  example, 
where  weight  and  volume  are  valuable  quantities  which  can  be 
directly  equated  to  dollar  cost.   An  unnecessarily  large  sup- 
port component  impacts  either  payload  or  endurance. 

It  is  certainly  conceivable  that  a  heat  exchanger  could 
be  designed  with  (a)  the  least  acquisition  cost,  (b)  the  low- 
est projected  operating  expenses,  and  yet  (c)  have  the  high- 
est weight  or  volume  impact  on  the  parent  design  and  there- 
fore by  unacceptable  to  to  "impact  cost". 

The  estimation  of  dollar  cost  of  a  particular  heat  ex- 
changer design  is  difficult  to  estimate,  primarily  due  to  the 
proprietary  nature  of  the  required  cost  estimating  relations. 
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Companies  that  design  and  manufacture  heat  exchangers  are  of- 
ten quite  reluctant  to  release  this  data. 

A  method  of  dollar  cost  comparison  will  be  proposed,  but 
due  to  lack  of  data,  will  not  be  fully  investigated. 

The  term  "cost",  therefore,  as  used  in  this  paper,  will 
most  generally  refer  to  "impact  cost"  and  will  serve  as  a 
means  of  relative  ranking  of  proposed  designs  specifically  by 
required  volume,  pumping  power  and  NTU  comparison. 

C.  Surface  Testing  and  Data  Presentation 

In  order  to  provide  the  required  friction  and  heat  trans- 
fer data  for  a  specified  finned  tubular  surface,  experimental 
testing  is  required. 

The  experimental  procedure  employed  by  Kays  and  London 
is  clearly  outlined  in  reference  [2], 

Regardless  of  the  technique,  the  results  include  fanning 
friction  factor  f,  Prandtl  number,  Reynolds  number,  and  Col- 
burn  j -factor.   Figure  2  is  an  example  showing  the  tradition- 
al method  of  data  presentation. 

D.  The  Task  of  the  Designer 

Assuming  that  the  designer  of  a  heat  exchanger  has  de- 
fined such  constraints  as: 

(a)  allowable  pressure  loss, 

(b)  temperature  change, 

(c)  amount  of  heat  to  be  transferred, 
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FIGURE   2:    Sample  data  plot   showing   traditional  method  of   data 
Presentation 
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(d)  volume  and  weight  restrictions, 

(e)  fouling  and  corrosion  considerations 

(f)  materials 

he  is  now  ready  to  select  the  heat  transfer  surface (s)  to  be 
employed. 

Having  compiled  a  great  number  of  candidate  surfaces 
from  (a)  existing  data,  or  (b)  having  designed  and  tested  his 
own  surface,  he  is  now  faced  with  the  task  of  comparing  these 
surfaces  to  determine  which  should  be  selected  for  the  specif- 
ic application  under  consideration. 

Figure  3  presents  two  examples  from  Ref  [2]  showing  the 
traditional  method  of  presenting  f  and  j  verses  RE  data  .  Note 
that  in  each  case,  Surface  B  has  a  superior  heat  transfer 
characteristic  at  a  given  Reynolds  number;  however,  Surface  A 
has  a  lower  (and  therefore,  superior)  friction  factor  at  that 
same  Reynolds  number.   Which  surface  is  "optimum"? 

There  appears  to  be  no  obvious  answer  to  this  question 
and  this  fact  has  prompted  several  methods  of  surface  compar- 
ison to  be  developed. 

It  is  intended  that  the  method  presented  herein  will  pro- 
vide a  solution  to  this  dilemma. 
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II.  PROPOSED  COMPARISON  TECHNIQUE 

A.  Derivation  of  Basis  of  Comparison 

Soland,  in  reference  [5] ,  provides  a  detailed  derivation 
of  the  proposed  surface  comparison  technique  used  to  analyze 
the  performance  of  crossflow  plate-finned  heat  exchangers. 

Theory  and  principle  inherent  with  the  Soland  method  are 
used  to  derive  a  comparison  technique  to  analyze  the  perfor- 
mance of  crossflow  finned  tubular  heat  exchangers. 

Comparison  of  the  performance  of  various  finned  tubular 
surfaces  assumes  the  following  quantities  are  held  constant: 

A)  W,  flow  rate 

B)  T,    ,  hot  fluid  inlet  temperature 

n,  in 

C)  T   .  ,  cold  gas  inlet  temperature 

Other  important  assumptions  are: 

E)  The  heat  transfer  resistance  of  the  tube 
walls  is  negligible. 

F)  The  controlling  heat  transfer  resistance 
is  assumed  to  be  on  the  finned  side  of  the 
heat  exchanger. 

Traditionally,  data  is  presented  in  terms  of  heat  trans- 
fer Colburn  modulus,  j,  and  friction  factor,  f,  based  on  the 
total  exposed  area,  A_,  as  a  function  of  Reynolds  number,  RE, 
based  on  minimum  free  flow  area,  A  ,  and  a  hydraullic  diameter, 

C 
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D. ,  of  the  flow  passage. 

The  proposed  comparison  method  converts  these  j  and  f  mag- 
nitudes to  new  quantities  jN  and  f„  based  on  the  heat  transfer 
area  of  the  bare  base  tube  surface,  A,  ,  and  a  Reynolds  number, 
RE  ,  which  is  referenced  to  the  open  flow  passage  area,  A  ,  as 
though  the  fins  were  not  present.   The  effect  of  the  fins  is 
accounted  for  as  an  increased  heat  flux  and  hence  larger  h  for 
for  the  bare  base  tube  surface  area.   In  order  to  incorporate 
the  effect  of  the  fins  into  j„,  the  metal  conductivity  must  be 
specified. 

Table  I  shows  the  proposed  new  definitions  and  compares 
them  with  the  analogous  definition  used  by  Kays  and  London  [2] 
Two  cases  must  be  considered  as  illustrated  in  Figure  4 . 

In  order  to  convert  data  presented  in  the  format  of  ref- 
erence [2]  to  the  new  basis,  various  ratios  are  derived  from 
basic  definitions,  equations  (1)  through  (10),  and  Figure  4. 
These  ratios  are  presented  in  Table  II. 

Two  assumptions  and  an  explanation  of  the  calculation  of 
radial  fin  efficiency  are  required  in  order  to  solve  for  the 
proper  fin  efficiency,  n,,  needed  in  the  conversion  relations: 

1.  The  results  presented  herein  will  assume  the  gas 
flowing  over  the  finned  tube  banks  is  air  whose 
thermodynamic  properties  will  be  evaluated  at  90 °F. 

2.  Some  of  the  results  presented  herein  specify  copper 
fins  [K*222  BTU/HR-FT-°F]  and  others  specify  Alumin- 
um fins  [K-117.6  BTU/HR-FT-°F] . 
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CASE  I;  Minimum  Free  Flow  Area  Occurs 
In  Spaces  Transverse  To  Flow 


♦Minimum  Free  Flow 
Dimension 


FLOW 


AF   =   atLN 


A.     =   ttD  LN  G„   =   W/A_   =   W/a. LN 

DO  N  r  t 


VT   =    FbLN 


4   AFb        4    atb        4    DoXL[XT-l] 
DN   =  —A—" 


7T  D 


7T 


XLD0  * 


CASE  II:  Minimum  Free  Flow  Area  Occurs 
In  Diagonal  Spaces 


♦Minimum  Free  Flow 
Dimension 


A_,=  2  a,LN 
F      d 


t=> 


A.  =  ttD  LN 
b     o 


VT  -  FbLN 


4  A„b 
DN  =  -A^ 


FLOW 


8  adb  _  8  "oW1' 

ttD  tt 


FIGURE  4 :  The  Two  Cases  of  Minimum  Free  Flow  Area  Occurance 

and  the  Associated  Geometric  Relations  of  Each  Case 
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3.  Figure  5  explains  a  proposed  generalized  procedure  for 
estimating  radial  fin  efficiency.   This  procedure 
provides  a  simplifed  method  of  efficiency  calculation 
in  lieu  of  using  more  technically  correct  Bessel  func- 
tion  type  relations. 
Figure  6  shows  an  example  of  data  presented  on  both  basis. 

It  should  be  noted  that  additional  curves  of  j..  verses  RE.. 

N  N 

would  result  for  other  magnitudes  of  fin  thermal  conductivity. 
For  any  heat  exchanger,  the  fluid  pumping  power  per  unit 
volume  on  one  side  is  given  by: 


Vos"  0   Vos 


From  equation  4b: 

AP,  =  £N  *    N  b  GN 

From  equation  2b: 


W  -  GN  Ap 


V         2     4 
os  g^  p    D„  v^e 

O         N    OS 


(17) 


(18) 
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0 


Compute  fin  efficiency  assuming  a 
Longitudinal  fin  of  Rectangular 
Profile i 

Tanh  ml 


where : 


2  h. 


jnr 


.50 


what  geometric  shape  describes  the 
fin  being  employed  on  the  tube 
surface? 


Radial  Pin  of  Rectangular 
Profile 


Tapered  (Hyperbolic)  Rad- 
ial fin 


an 


correction 


CA  -  DA  lnfn^) 


where » 

r. 
CA«  -.00259  -  .0084  in  -°- 


D.-  -.02426  -  -4521  In  -^ 


(Ref  3  page  264) 


Rectangular  Profile  fin 

sheet 
[see   below  for  I   defn] 


Correction" 


[Ref  4  page  3-116) 


An 


correction 


Same  as  Radial  Fin  of 
Rectangular  Profile 
using  rg( staggered  bank) 

50   XdDo 
d  o 


-:  -  m 


[Ref  4  page  3-116] 


71  i  -corrected 


n«-t 


-  An 


correction 


Compute  Surface  ns 
».  '  l  '   *^U  -  ^-corrected' 


FIGURE   5;      Computing   fin  efficiency  of  a  Radial   fin 
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FIGURE  6;  Comparison  of  Colburn  j  factors  (j  and  j  ) ,  and  com- 
parison of  friction  factors  (f  and  f  )  for  the  Sur- 
face of  Figure  94  [Ref  2] . 
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For  the  same  fluid  at  the  same  temperature  level,  y  and  p 
are  constant: 


P      fN  REN  AF  b 

P_   «       N  (19) 

os      D*  v 

N    OS 


Figure  7  shows  the  geometry  involved  in  the  calculation 
of  [A  b]/V   for  both  case  I  and  case  II. 


The  heat  transfer  for  any  heat  exchanger  is  given  by: 


S  =  elTh  in  "  T^  iJw  c«  (20) 

n ,  in    c  $  in    p 


For  any  given  flow  arrangement,  e  can  be  related  to  NTU 
by  either  a  closed  form  equation  or  graphical  curve  similar  to 
that  shown  in  figure  8,  where 


NTU  H  ^  (21 


The  relationship  between  e  and  NTU  is  always  monotonical- 
ly  increasing.   In  other  words,  an  increase  in  A  h„  causes  an 
increase  in  NTU  which  means  e  and  hence  q  are  greater,  given 
that  the  fluid  properties  and  flow  rate  are  held  constant. 

From  equations  (21)  and  (7b) : 
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Case  I :  Minimum  Free  flow  area  Transverse  to  Flow  Direction 


Total  Volume  on  one  Side: 


Vog=  [XTD0HXLD0]L   -  JD?  L 


L  =  dimension  into  paper  along 
flow  length  of  tubes. 


"Ap  b"  Volume:  Y////A 

A_  b=  [X_D   -  D  ]L  XT  D 
r      T  o    o    L  c 


Ratio: 


V. 


OS 


[xT  -  i]xL 

XTXL  -  it/4 


Case  II:   Minimum  Free  Flow  area  in  the  Diagonals 


Total  Volume  on  one  Side: 
Vos  ■  [Do  XD  8ina  -  J°o' 


'*F 


b"   Volume: 


A,   b   =    2    XL   DQ    [XD   -    1] 


Ratio: 


AF   b        4[Xp   -   1]XL 

os       X_sina  -  x 

EL 1— 


v^ 


9 


M>" 


FIGURE  7:   Calculation  of    tAF  bl//Vos  for  both  case   I  and  II. 
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NTU  .  Ab  hN    .  *b  jN  GN  Cp 

V     V   w  c    ■>  * 
os    os    p  P2  3  V   W  C 
r   r   os    p 


From  equations  (3b)  and  (lb) 


NTU  =  4  AF  b  JN  ^N  y 

VOS    D?,   P2  3  W  V 

N  r     os 


(22) 


(23) 


or 


AhN 

4    C       M 
P 

JN    REN  AF   b 

Vos 

p2  3 

r 

2          V 

(24) 


For  the  same  fluid  at  the  same  temperature  level,  C  ,  y , 
and  P  are  constant: 


A  hN  .   ^  ^N   AF  b 

v —   — —    

N       os 


And  for  the  flow  rate,  W,  also  constant: 


NTUa   jN  ^N   AF  b 
V  2 

N      OS 


(25) 


(26) 
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NTU 


FIGURE  8:   Typical  Graphical  Relationship  between  heat  ex- 
changer effectiveness  and  NTU.   In  general  it  is 
possible  to  express: 

e  =  *[NTU,  C  .  /C  m.    flow  arrangement] 
min  max 
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Recognizing  that  the  f  and  j  verses  RE  plots  for  flow 
over  finned  tube  banks  from  reference  [2]  can  be  closely  ap- 
proximated by  a  straight  line  on  a  log-log  plot,  implies  that 
f  and  j  can  be  represented  by  an  equation  of  the  form: 

Bf 
f (•)  -  Af  «■ 

B. 

j,  x  =  A.  RE   D 
J  (e)     J 

Table  I  of  Appendix  III  lists  the  linearized  equations  of 
f  and  j  for  the  surfaces  represented  by  figures  92  -101  of 
reference  [2] .   This  simplification  allows  for  ease  in  pro- 
gramming the  derived  surface  comparison  equations  as  demon- 
strated in  the  procedure  outlined  in  Appendix  I. 

In  any  event,  the  performance  parameters  may  be  calculated 

from  enhanced  surface  data  once  the  data  is  presented  in  the 

form  f„  verses  RE„  and  j„  verses  REXT.   This  presentation  may 
N  N      JH  N 

be  (a)  graphical  due  to  the  complex  relations  presented  by  f 

and  j  verses  RE,  as  suggested  by  Soland  or  (b)  analytical,  as 

outlined  in  Appendix  I,  due  in  large  part  to  the  linear  nature 

of  f  and  j  verses  RE  on  a  log-log  plot. 

B.  Method  of  Surface  Comparison 

Having  expressed  the  enhanced  surface  data  in  the  form  fN 

verses  RE„  and  j„  verses  RE„,  it  is  a  relatively  simple  matter 
N       N  N 
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to  calculate  the  performance  parameters  of  equations  (19)  and 
(26)  and  plot  them.   Figure  9  is  an  example  of  such  a  plot. 
The  performance  parameters  of  two  hypothetical  surfaces,  label- 
ed 1  and  2,    have  been  plotted  in  order  to  demonstrate  how  a 
determination  of  heat  exchanger  relative  performance  may  be 
made. 

Four  different  comparisons  are  immediately  available  from 
figure  9  and  are  indicated  by  points  a,  b,  c,  and  d  on  surface 
2.   Point  o  on  surface  1  represents  the  reference  heat  exchanger 

design  to  which  each  of  the  four  points  on  surface  2  will  be 
compared . 


Case  A;   Same  heat  exchanger  shape  and  volume  (L  =L  ,  V  =V  , 

A_  =A_  ) .   Because  W  and  A_  are  fixed : 
a   o 

lDN  ] 
RE„  -  RE„      a  (27) 

a      °  [DN  ] 
o 

The  results  of  this  comparison  are  easily  obtained  as 

the  ratios  of  ordinate  values  [NTU/V  1  and  abscissa 

values  P/V   of  the  two  surfaces. 
os 


Case  B;   Same  heat  exchanger  volume  and  pumping  power  (V]3=v0' 

P.  =P  ) .  Point  b  is  located  on  a  vertical  line  throug 
o  o 

point  o  since  the  pumping  power  per  unit  volume  is 
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FIGURE  9;  Performance  Parameter  Curves  For  Two 
Hypothetical  Surfaces  Showing  Points 
Used  in  Sample  Comparisons. 
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Case  B;  (continued) 

identical  in  both  exchangers.   The  NTU  ratio  of  the 

two  heat  exchangers  is  obtained  as  the  ratio  of  the 

ordinate  values  from  Figure  9. 

Case  C;   Same  pumping  power  and  number  of  transfer  units. 

(P  =P  ,  NTU  =NTU  )  Point  c  is  located  on  a  line  hav- 

C   O       CO 

ing  slope  equal  to  unity  and  through  point  o  since 
both  NTU  and  P  are  constant  and  each  axis  is  inverse- 
ly proportional  to  volume.   The  ratio  of  the  volume 
required  using  surface  2  to  the  volume  required  using 
surface  1  is  the  ratio  of  either  ordinates  or  abscis- 
sas at  points  o  and  c,  respectively. 

Since  material  costs  are  usually  a  function  of  volume 
or  material  weight,  and  thus  heat  transfer  surface 
area,  it  is  through  Case  C  that  a  relative  cost  com- 
parison may  be  attempted: 

[Total  Cost]2    [$/FT2]2  [Bl2  [VT] 2 


[Total  Cost^    [$/FT2]1  [$]x  [VT11 


Note  that  as  long  as  surface  2  lies  above  surface  1, 
the  result  will  be  a  smaller  volume  required  to  do 
the  same  "job",  [see  Figure  9] 


Case  D:   Same  volume  and  number  of  transfer  units.   (Vd=Vo' 

NTU j^NTU  or  q  =q  ) .   This  case  compares  the  pumping 

uO      a   O 
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Case  D;  (continued) 

power  required  by  surface  2  to  that  required  by  sur- 
face 1  for  the  case  when  both  heat  exchangers  yield 
the  same  overall  heat  transfer  performance.   The  com- 
parison is  made  by  taking  the  ratio  of  the  abscissas 
(points  o  and  d)  in  Figure  9,  since  a  horizontal  line 
on  the  performance  plot  has  a  constant  value  of 

NTU/V  .   The  surface  2  Reynolds  number  for  this  case 
os 

is  the  smallest  of  all  four  cases  and  consequently, 
the  flow  area  is  the  largest.   Since  the  volume  is  the 
same  as  the  surface  1  exchanger,  the  length  must  de- 
crease. 


C.  Summary 

Using  the  equations  derived  herein  it  is  possible  to  con- 
struct a  performance  parameter  plot  of  NTU/V   verses  P/v 

OS  OS 

(see  Figure  9) ,  from  which  it  is  possible  to  obtain  useful 
performance  comparisons  between  two  heat  exchangers  for  four 
different  criteria. 
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III.  COMPARISON  OF  FINNED  TUBE  SURFACES 

A.  General 

Using  the  procedure  derived  in  Chapter  II  and  outlined  in 
Appendix  I,  all  the  finned  tube  surfaces  were  plotted  as  per- 
formance curves  of: 


*H  ^N  *F  b      _      fN  MN  AF  b 
x— verses   ■? 

DNVos  DN  Vos 


These  performance  curves  assume  constant  W,  T   .  ,  and 

c ,  in 

Tu  «_#  and  are  equivalent  to  NTU/V  „  verses  P/V 
n,in  ^  os  os. 

Figure  10  shows  the  resulting  performance  parameter  plot 
presenting  only  the  highest  performance  curve  at  each  nominal 
diameter  [DN1 .   This  study  considered  sixteen  different  sur- 
faces having  eleven  different  nominal  diameters.   Calculated 
performance  results  are  presented  for  each  surface  considered 
in  Appendix  II. 

B.   Comparison  to  Smooth  Surface 

All  of  the  sixteen  finned  tubular  surfaces  were  then  com- 
pared to  a  smooth  surface  having  no  enhancements  (fins)  and  a 
nominal  diameter  identical  to  that  of  the  enhanced  surface. 
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PERFORMANCE  PARAMETER  CURVE  FOR  FINNED  TUBULAR  SURFACES 


FIGURE  10 
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Note  that  subscript  e  will  be  used  to  denote  the  enhanced  sur- 
face and  the  subscript  s  will  be  used  to  denote  the  smooth 


surface,    i.    « 

1-1 


using 


&e|zi8ik< 


■K 


Snjpogh  surface  -*s   N  HF   calculations  were  completed 
0  e  w 

~  1 1 


DN  Vos 


turbul«itL.f  Jbw  o  rer  staggered  tube* 


au  ;kas  [7]  correlation  for  forced-convection, 


banks  with  the  appropriate 


nominal;  diameter 

t   S   - 

8  :  | 

"  Js%     m    i 

!»      *>     %< 
D     A      ft 

1$  ord£r  to 

*     I     -     . 


predict  sjhoSth" ttbe  bank  ColbUfn  Modulus    [j    ]  ,    equation    (28) 


*    *• 


I*.         w 


suits  SireJ  shpvlh 


Note  tfcatr;  j 

ced  surface,  in 


[DN1. 


[RE]"0'396  [XT/XL] 


/6  [Do/DN] 


-0.400 


(28) 


justify  the  use  o4  the  above  relation  to 


• 


was  used  ttocpx&d  Let  the  actual  Heat  transfer  data   for   flow 
*     •     I     «  • 

over  bare£  ttfbet  binks  as  presented  i  it  reference    [2].    The  re- 


•  j»  * 


m  figure  11.   Maximum  deviation  between  the 


experimental  valie  of  t   (from  reference  2)  and  the  calculated 
value  of  »j  ''from 


equation  (28)  is  ^7.5%. 

may  be  presented]  for  each  individual  enhan- 
the  ffcrm:   *      *      • 
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(29) 
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N        W  N 
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Equation  (26)  may  then  be  written  for  the  enhanced  sur- 
face as  "if  smooth": 


NTOS_    Us}     RENAFb  .,„. 

V  2  K       ' 

OS        D„    V 

N     OS 


3 

f   RE   A   b 

Smooth  surface   s   N  F  '   were  completed  using  a  corre- 

DN   Vos 
lation  for  friction  factor  (f  )  derived  from  first  principles 
using  the  Kays  and  London  [2]  data  for  turbulent  flow  over 
staggered  bare  tube  banks  as  the  data  base: 


Case  I: 

f   =  0.7184  RE"0'2407  [X^-l]0-703  IX,]0*             (31) 

S  1               Jj 

where:  CT  =  0.828  -  16.596  RE*0,4277 


Case  II; 

C 

f   =  1.0000  RE"0'1892  [X  -I]0*703  [XT]  D  (32) 

S  O  Jj 

where:  CQ  -  130.10  RE-0,7284 

Reference  to  Case  I  and  Case  II  differentiates  between 
minimum  free  flow  area  occur ing  transversely  (I)  or  diagonally 
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(II)  to  the  flow  direction,  (see  Figure  4) . 

In  order  to  justify  the  use  of  equations  (31)  and  (32)  to 
predict  smooth  tube  bank  friction  factor  (f  ) ,  these  relations 
were  used  to  predict  the  actual  experimental  friction  factor 
data  for  flow  over  bare  tube  banks  as  presented  in  reference 
[2J .   The  results  are  shown  in  Figure  12.   Maximum  deviation 
between  the  experimental  value  of  f   (from  reference  2)  and 
the  calculated  value  from  either  (31)  or  (32)  is  =3.0%. 

The  friction  factor,  for  each  individual  enhanced  surface 
may  be  presented  in  the  form: 


f   =  A      REBf(8) 
s    Af(s)  ^N  (33) 


Equation  (19)  may  then  be  written  for  the  enhanced  sur- 
face as  "if  smooth": 


f   RE.*  A   b 
f-=   S   4  (34) 

DN   Vos 


Of  the  sixteen  surfaces  investigated,  five  have  a  nominal 
diameter  [D„]=1.60  inches.   Since  this  offers  an  excellent 
means  of  comparison,  cases  A,B,C,  and  D  where  calculated  for 
these  five  surfaces: 
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Figure  13;   Shows  the  ratio  of  ordinates  [NTU  /NTU  ] 

verses  RE.,  for  Case  A.  [Same  heat  exchanger 
shape  and  volume]   Since  both  the  smooth 
and  enhanced  surfaces  have  identical  nomin- 
al diameters,  the  nominal  Reynolds  number 
is  also  identical  for  each  surface,  [see 
equation  27;  results  read  from  Appendix  II] 


Figure  14 :   Shows  the  ratio  of  abscissas  [P  /P  ]  verses 

RE.,  for  Case  A.  Since  both  the  smooth  and 
N     

enhanced  surfaces  have  identical  nominal 
diameters,  the  nominal  Reynolds  number  is 
also  identical  for  each  surface.   [see 
equation  27; results  read  from  Appendix  II] 


Figure  15:   Shows  the  ratio  of  ordinates  [NTU  /NTU  1 

■  e     s 

verses  RE„  for  Case  B.  [P=constant  and 

N       

V=constant] 


Figure  16:   Shows  the  ratio  of  abscissas  or  ordinates 

to  produce  [V  /V  ]  verses  RE.,.   Recall  that 
each  axis  of  Figure  10  is  inversely  pro- 
portional to  volume.  This  is  Case  C. 
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Figure  17:   Shows  the  ratio  of  abscissas  [P  /P  ] 

verses  RE.,  for  Case  D.  [NTU=constant,  and 

N       

V^constant] 

Higher  ratios  are  preferred  for  Figures  13,15,16,  and  17 
while  a  lower  ratio  is  most  desirable  in  Figure  14. 

As  a  matter  of  secondary  interest,  of  the  five  surfaces 
with  identical  nominal  diameters,  there  are  three  discrete  val- 
ues of  the  ratio  [A_  b]/V   : 

F      OS 

Surface  [AF  b]/Vos 
92  &  93  1.4282 

96  &  97A  1.1859 

100  1.4014 

Comparison  of  many  surfaces  having  identical  nominal  dia- 
meters to  a  single  common  smooth  surface  having  the  same  nom- 
inal diameter  permits  the  best  relative  comparison  between  en- 
hanced surfaces. 
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FIGURE  17;  Performance  comparison  results  for  CASE  D  [i.e.  NTU 
and  V  are  constant]  D„  -   1.6 
enhanced  and  plain  surfaces. 


and  V  are  constant]  D„  -   1.60  inches  for  both  the 
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IV.  HEAT  EXCHANGER  DESIGN  PROBLEM 

A.  Description 

Appendix  III  is  a  procedure  which  can  be  used  for  sizing 
crossflow  finned  tubular  heat  exchangers.   The  following  data, 
taken  from  reference  [2] ,    is  used  to  determine  the  required 
heat  exchanger  size  for  the  given  conditions. 

The  design  is  that  of  an  intercooler  for  a  5000  SHP  gas 
turbine  plant. 


DESIGN  CONSTRAINT 

SIDE  I 

SIDE  II 

Surface 

Smooth  Tube 

Finned  Tube  [98C] 

Dh 

0.0625  FT 

0.0445  FT 

6 

N/A 

0.0010  FT 

S 

5.9291  FT2/FT3 

61.9  FT2/FT3 

VAT 

N/A 

0.8350 

I 

N/A 

0.3445  IN 

W 

400,000  LB/HR 

200,667  LB/HR 

T. 
in 

60°  F 

260°  F 

T 
out 

82°  F 

80°  F 

AP 

1.554  PSI 

0.3060  PSI 

P 

2.360  LB/HR-FT 

0.0482  LB/HR-FT 

S 

1.000  BTU/LB-°F 

0.2436  BTU/LB-°F 

P« 

62.38  LB/FT3 

0.1755  LB/FT3 

m 
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DESIGN  CONSTRAINT 

SIDE  I 
6.8000 

SIDE  II 

P 

r 

0.7000 

K 
m 

N/A 

222  BTU/HR-FT-'F 

F 

2.725  IN 

N/A 

b 

1.750  IN 

N/A 

A 

0.4418  IN2 

Calculate 

Using  the  design  procedure  outlined  in  Appendix  III,  the 
principle  dimensions  are: 

X  =  6.4  3  FT 
Y  =  26.11  FT 
Z  =  0.63  FT 
Total  Volume  =  105.46  FT3 

This  is  a  possible  heat  exchanger  design  that  satifies 
the  given  constraints. 

B.  Attempt  to  Minimize  Total  Volume 

If  minimization  of  total  volume  is  a  required  design  con- 
straint, the  engineer  has  several  options  the  limits  of  which 
might  be  bounded  by: 

(A)  Design  a  heat  exchanger  to  the  given  specifications 
using  each  known  surface  for  which  f  and  j  data  are 
available.  In  this  study  sixteen  surfaces  are  con- 
sidered. 

(B)  Noting  the  marked  superiority  of  reference  [2]  figure 
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(B)  continued 

94,  apply  the  method  developed  in  Chapter  II  to  es- 
timate total  volume  savings  relative  to  the  base  or 
reference  heat  exchanger  designed  in  part  I. 

Option  (A)  will  be  eliminated  as  being  impractical  due  to 
the  required  effort. 

Option  (B)  will  be  shown  to  be  the  most  direct  and  highly 
accurate  approach. 

The  f  and  j  verses  RE  data  for  reference  [2]  figure  98C 
and  94  are  presented  in  Figure  18.   The  problem  presented  cor- 
responds to  Case  C  described  in  Figure  16. 

Using  the  heat  exchanger  design  procedure  outlined  in 
Appendix  III  and  replacing  the  finned  tube  side  with  the  sur- 
face of  figure  94  [Ref  2]: 

X  =  1.50  FT 
Y  =  10.89  FT 
Z  ■  1.76  FT 
Total  Volume  =  28.68  FT3 
A  73%  reduction  in  volume  would  result  by  employing  the 
surface  of  figure  94  [Ref  2] . 

In  order  to  use  the  comparison  method  developed  in  Chap- 
ter II  to  quantitatively  predict  the  volume  savings,  a  perfor 
mance  parameter  plot  is  constructed  for  each  surface  consider- 
ed.  Figure  19  presents  the  only  two  surfaces  being  considered 
in  this  example. 
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FIGURE  18:  f  and  j  verses  RE  data  for  the  two  surfaces  con- 
sidered in  the  design  example. 
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As  outlined  in  Chapter  II  for  Case  C: 

(1)  From  the  design  point  of  the  reference  heat 
exchanger  (Figure  98C)  graphically  or  analyt- 
ically construct  a  line  of  slope  1. 

(2)  The  intersection  of  this  line  with  the  new 
surface  (Figure  94)  will  yield  the  performance 
parameters  at  the  design  point  of  the  new 
heat  exchanger. 

(3)  The  predicted  volume  reduction  is  found  from 
the  ratio  of  either  the  ordinates  or  abscissas 
of  the  two  surfaces. 

The  performance  parameters  of  the  reference  heat  exchang- 
er design  are:  [Reference  2  figure  98C] 


^  =  401.20  [1/IN2] 

OS 


P  12        4 

=  2.93  X  10x*  [1/IN  ] 


V 
OS 


The  performance  parameters  of  the  heat  exchanger  that 
would  result  from  using  the  surface  described  by  reference [2] 
figure  94 : 

^  =  1634.68  [1/IN2] 
os 


P  13        4 

=  1.19  X  10XJ  [1/IN4] 
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In  other  words,  a  predicted  volume  reduction  of  75%  is 
anticipated  using  the  comparison  method  developed  herein. 

C.  Summary 

Two  surfaces  have  been  compared  from  a  minimization  of 
total  volume  point-of-view.   It  has  been  shown  that  the  per- 
formance parameters  developed  in  Chapter  II  allow  us  to  com- 
pare any  number  of  surfaces  without  going  through  a  complete 
heat  exchanger  design  calculation.   Where  there  may  well  be 
hundreds  of  different  surfaces  to  be  considered,  the  value  of 
such  a  comparison  technique  cannot  be  overestimated. 
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V.   CONCLUSIONS 

1.  The  comparison  method  of  Soland  [5]  was  used  as  a  basis 
in  developing  a  comparison  method  for  finned  tube  heat  ex- 
changer surfaces. 

2.  Heat  Exchanger  Performance  can  be  compared  on  four 
different  bases: 

a.  Same  shape  and  volume  of  heat  exchanger 

b.  Same  exchanger  volume  and  pumping  power 

c.  Same  pumping  power  and  NTU 

d.  Same  volume  and  NTU 

3.  All  of  Kays  and  London  [2]  finned  tube  surfaces  and 
several  surfaces  supplied  by  the  Trane  Corporation  were  com- 
pared on  the  above  basis. 

4.  The  method  developed  herein  was  applied  to  a  practical 
heat  exchanger  design  problem  and  found  to  accurately  predict 
relative  volume  savings  between  different  surfaces. 
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APPENDIX  I 

A.  Purpose 

A  frequent  criticism  of  academic  solutions  to 
"optimization"  of  design  is:   "So  what,  .  .  .  how  does  this 
help  the  engineer  in  industry?"  [8] 

It  is  hoped  that  this  section  will  give  the  "engineer 
in  industry"  the  needed  background  and  technique  to  effective- 
ly compare  finned  tubular  heat  exchanger  surfaces. 

A  step-by-step  flow  diagram  using  the  performance 
comparison  equations  derived  in  the  text  of  the  thesis  is 
presented  in  Figure  A-l.   This  procedure  is  easily  programmed 
for  use  on  either  a  (1)  standard  computer  or  (2)  programmable 
calculator. 

In  keeping  with  the  spirit  of  the  criticism  of 
Reference  [8] ,  a  program  for  use  on  a  personal  programmable 
calculator  will  be  presented  herein.   It  is  intended  that 
this  simplified,  easy-to-use  approach  will  demonstrate  that 
the  "engineer  in  industry"  needn't  resort  to  often  costly 
computer  time/programs  in  order  to  effectively  "optimize" 
a  design. 

B.  Proposed  Program(s)  for  Calculating  Surface  Performance 
Parameters  on  a  Programmable*  Calculator 

The  following  programs  follow  the  general  outline 
of  Figure  A-l: 

Figure  A-II:   For  calculating  Surface  Performance 
Parameters  for  the  case  in  which  the 
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CALCULATION  OF  FINNED  TUBULAR  HEAT  EXCHANGER  SURFACE  PERFORMANCE  PARAMETERS 
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FIGURE  A-l:    Performance  Parameters 
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minimum  free  flow  area  occurs  transver- 
sely to  the  flow  direction. 

Figure  A-II-A;   Definition  of  geometric  constraints/ 

constants  to  be  stored  in  memory. 

Figure  A-II-B:   Sample  output  for  transverse  case. 


Figure  A-III;    For  calculating  Surface  Performance 

Parameters  for  the  case  in  which  the 
minimum  free  flow  area  occurs  in  the 
diagonals.   [see  Figure  4] 

Figure  A-III-A:  Definition  of  geometric  constraints/ 

constants  to  be  stored  in  memory. 

Figure  A-III-B:  Sample  output  for  diagonal  case. 

C.  Operation 

Step  Description 

1.  Select  and  Load  Program 

2.  Store  Constants  as  outlined  in  either 

Figure  A-II-A  or  A-III-A. 

3.  Press  Key  A 

4.  Compile  Output 

5.  Stop  Program  when  range  of  RE  data  is 

exceeded . 
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PROGRAM*  FOR  CALCULATING  SURFACE  PERFORMANCE  PARAMETERS 
FOR  A  CROSSFLOW  FINNED  TUBULAR  HEAT  EXCHANGER  IN  WHICH 
THE  MINIMUM  FREE  FLOW  AREA  OCCURS  TRANSVERSLY  TO  THE  FLOW 

DIRECTION 
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FIGURE      A-II-A 
GEOMETRIC   CONSTRAINTS/CONSTANTS    TO   BE    STORED    IN   MEMORY 

[Transverse  Case] 
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FIGURE  A-II-B 
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FIGURE  A-III 


PROGRAM*  FOR  CALCULATING  SURFACE  PERFORMANCE  PARAMETERS  FOR  A 
CROSSFLOW  FINNED  TUBULAR  HEAT  EXCHANGER  IN  WHICH  THE  MINIMUM 
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FIGURE   A-III-A 


GEOMETRIC    CONSTRAINTS/CONSTANTS    TO   BE    STORED    IN   MEMORY 

[Diagonal   Case] 


no 

or* 

DESCRIPTION 

MO 

0T» 

DESCRIPTION 

• 

«D 

■pacifiad  diagonal  tub*  pitch. 
Ratio  of  diagonal  tuba  (pacing 
to  th*  tuba  outalda  diameter. 

10 

V*T 

Ratio  of  fin  area  to  total 
area. 

X 

Af«*> 

from  linearized  aquation  for 
friction  factor  (f)  of  the  en- 
hanced aurfaca  of  the  font 
A  RE  * 

11 

t//Z 

Ratio  of  fin  length  to  fin 
thlcknaea.  U//Inch| 

a 

'«<•> 

rroa  linearized  aquation  for 
friction  factor  (f)  of  tho  en- 
hanced aurfaca  of  the  formi 
A  RE  ' 

12 

i 

Property  conatant  defined  byi 
(Cp  III/  (ltB  t'/'l          (ND] 

'  i 

*«(•> 

rroa  linearized  aquation  for 
friction  factor  (f)  of  the  if 
■sooth  aurfaca  of  the  font 

an" 

13 

*!<•> 

froa  linoarized  equetlon  for 
Colburn  modulus  (j)  of  the  en- 
hanced aurfaca  of  tha  forai 
A  RE  B 

4 

■«■> 

Froa  linaarltad  aquation  for 
friction  factor  (f)  of  the  if 
•Booth  aurfaca  of  the  fora: 
A  AE» 

14 

Bl<*> 

From  linearized  equation  for 
Colburn  modulus  (j)  of  the  en- 
hanced aurface  of  tha  form: 

Alt' 

1 

D2h 

•pacifiad  hydraullic  dlaaatar 
of  the  flow  paaaaga  Ilnchaa) 

IS 

«T 

of  tranavaree  tuba  apacing  to 
tuba  outalda  diameter. 

• 

'x 
*a 
*j 
*4 

Conatanta  derived  froa  graph  on 
paga  2(4  of  Reference  (3| .   They 
•re  uaed  In  the  ooaputatlon  of  a 
correction  factor  applied  to 
h<_L  1»  ordar  to  permit  tha  uaa 
of  tha  alapller  Tanh  at  aquation 

for  fin  temperature  effeotlve- 
aoaa  In  lieu  of  aeaeel  function 
typo  aquatlona. 

14 

*)<a) 

rroa  linearized  equation  for 
Colburn  aodulua  (j)  of  tha  "If 
aaooth*  aurfaca  of  tha  forai 

are" 

7 

17 

»J(.) 

Froa  lineerized  equetlon  for 
Colburn  Hodulua  (j)  of  th*  ■if 
aaooth*  aurfaca  of  th*  forai 
A  RE» 

• 

IS 

Pin  ID 

Input  on* > 

0.000  If  fin  la  tapered 
-1.000  If  fin  ia  a  radial  fin 
of  rectangular,  profile 

t 

IS 

P 

Ratio  of  outalda  tub*  radlua 
to  outalda  fin  radlua. 

A 

Data  ID 

Any  daalrad  number  to  permit  tho 
ueer  to  Identify  tha  aurfaca  be- 
ing procaaaad. 

0 

Xl 

longitudinal  tuba  pitch.   Ratio 
of  longitudinal  tube  apacing  to 
tub*  outelde  diameter. 

■ 

1 

■pacified  heat  tranefer  area  par 
total  volume  ratio.  iln3/ln3| 

E 

a 

Angle  defined  by  longitudinal 
centerline  and  diagonal  cen- 
terllnei  (dcgreoel 

e 

D. 

•pacifiad  tuba  outalda  diameter. 
Ilnchaa) 

: 

' Increment 

Daalrad  Increment  of  RE  .   A 
value  of  10,000  ia  recoramneded . 
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FIGURE   A-IIIB 


SAMPLE  OUTPUT 
[diagonal   case] 


FIGURE  ID#< 


Vf 


[A_  b]/V  ■ 
F     os 


-^98.4000 

2.0413 

0.2173 
-0.2590 

1.1125 
-0.1971 

0.1904 

0.0027 
-0.0084 

0.2636 
-0.4521 

2.1214 

■^13.5882 

10.0860 

1.0336 

->1.5332 

991.4751 


0 
1 
2 
3 
4 
5 
6 
7 
8 
9 
A 
B 
C 
D 
E 
I 


Print  out  of  Regis- 
ters  0-9  and  A-I 


D 


N 


REN/RE 


Value  of  REM  being 


N 


evaluated 


0.8350 

3.1448 

0.0001 

0.0813 

-0.3264 

3.5207 

0.5965 

-0.3960 

-1.0000 

0.5290 


P  /V 
s   os 


NTUs/Vos' 


10000.0000 

3.743772534+10 

->1. 370584009+10 

320.2228 

>  52.9601 


20000.0000 

2.502842147+11 

9.564330240+10 

507.5132 

80.4954 


0 
1 
2 
3 
4 
5 
6 
7 
8 
9 


Print  out  of  Regis- 
ters  10-19 


P  /Vos 


-NTU  /V  a 
e   os 
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APPENDIX   II 


The  purpose  of  this  appendix  is  to  present: 

(1)  A  full  description  of  the  geometry  of  each  sur- 
face evaluated  in  this  study. 

(2)  Friction  factor  (f)  and  Colburn  Modulus  (j)  as 
a  funtion  of  Reynolds  number  (RE)  for  each  sur- 
face evaluated  in  this  study. 

(3)  A  complete  compilation  of  the  surface  perfor- 
mance parameters  required  to  compare  each  sur- 
face evaluated  in  this  study. 
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APPENDIX  III 
SIZING  CROSSFLOW  FINNED  TUBULAR  HEAT  EXCHANGERS  FOR  A  GIVEN  JOB 

Crossflow  finned  tubular  heat  exchangers  have  found  app- 
lication as  gas  turbine  plant  intercoolers,  air  conditioning 
unit  heat  exchangers,  aircraft  engine  coolers,  and  numerous 
other  uses.   These  applications  most  often  involve  gas-liquid 
service. 

The  "job"  of  the  heat  exchanger  will  be  defined  as  trans- 
fering  a  specified  amount  of  heat  between  two  fluids  at  given 
flow  rates  (W)  and  with  specified  amounts  of  pumping  power  on 
each  side.   The  following  values  are  specified: 

V  Tini'    Toutl'    AP1'    Pini'    Dlh'    Pr°PertieSl 

V  Tin2'   Tout2'   AP2'   Pin2'   D2h'   Properties., 

Subscript  1  will  refer  to  the  hotter  fluid  inside  the 
smooth  tubes  and  subscript  2  to  the  colder  fluid  flowing  over 
finned  tube  banks. 

While  core  pressure  drop  accounts  for  by  far  the  greatest 
percentage  of  the  total  pressure  drop,  entrance  and  exit  loss- 
es would  have  to  be  included  in  the  final  design.   [Ref  1  Chap- 
ter 4] 

Figure  A-l  shows  the  heat  exchanger  arrangement. 
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FIGURE  A-l 


1 


SIDE  1 


(1)  In  the  air-conditioning  industry, 
fins  which  are  common  to  all  tubes  in 
the  bundle  as  distinct  from  having 
each  tube  separately  finned  are  often 
used.  This  is  called  Plate  finning. 


U 


rnw/x    vmm       mm 


(2)  L-fins   (3)  Extruded   (4)  G-fins 
I 
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The  Basic  Equations 

The  first  portion  of  the  procedure  involves  determination 
of  A   ,  A  ,    X,  Yf  and  Z  based  on  the  specified  values  of  AP. 


cl    c2 


and  AP~  and  an  assumed  value  of  RE?.   Then  based  on  heat  trans- 
fer considerations,  it  is  determined  if  the  proper  heat  balance 
exists.   If  not,  another  value  of  RE-  is  assumed  and  the  process 
is  repeated. 


4  A. 


h   ' 


Dlh  Afr. 


-rr  D 


lh 


F  b 


for  round 
tubes 


D1U  =  D.    , 
lh    inside 


(1) 


4  A 


e„  = 


D2h  Afr. 


=  A  specified  geometric  constant  for 

surface  2;  See  Figure  A-II 

lD2h  "  4rh' 


(2) 


i-  L      4     i 

fr, 


(3) 


rr„      4 


(4) 


K     ■  t3i  — 1  lxzl  =  K,  X  z 

cx      1   4  1 


(5) 
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A^   =  13,  -=£]  IYZ]  =  K0  Y  Z  (6) 

C2      Z   4  ^ 


G  =  ^-  (7) 

A 
c 


G,    'V^1    [Ml  tAc2] 


G2    tW2/AC2]    [W2]  [ACil 


From  equations  (5)  and  (6) 


(8) 


G,    [W-]  [K-  YZ] 

_±  =  — i £ =  K  -  o: 

G2    [W2]  [K1  XZ1       Z 


W   K 
where  K^  = — = — - 
w2  Kl 


4p  =  4  f  Lg2  (10) 

Dlh  2  *o  PI» 


where  L  is  the  flow  length 
Side  1:  L  =  Y 
Side  2:  L  =  X 
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Gl  _  [AP1  Dlh  pml'  tf2]  [X1 
°!    liP2  D2h  pm2]  tfl]  [Y1 


n  <=■    .     1AP1  Dlh  Pml' 
Define  K„  =  -nr=; — =; — r 

4    [AP2  D2h  pm2] 


(11) 


From  Equations  (9)  and  (11) : 


or 


G?  tf  o3     [X]  2       m2 

4  =   K.    — - =   K*      ill  (12) 

g\      4  [fx]  [y]      J    txr 


f,         K*    [Y]3  m3 

-1  =   -J *.  =   K-    i^-  (13) 

f±         K4     [X]J  S     [X]J 

where  K,.   =  


Eliminate    [Y/X]    with  equations    (9)    and    (13) 


f3         [G.]3    [K_]  [G   ]3 

2    -         1         !_      =   k  - (I4) 

fl         tG2]3    lK3l3  6lG2]3 
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K5 
where  K,   =  —^ 

6   v3 
K3 


AP2  P2  Sl  W2 
AP1  pl  B2  Wl 


G  =  BO. 


Substitution  in  equation  (14)  yields 


—  =  K 


[RE1  "1  D2h  i 
6  [RE2  4   Dlh  ] 


(15) 


RE1  fl 


f2  ^2  ^2  Dlh 


K, 


3  ^3 
yl  D2h 


(16) 


Solving  for  RE, : 


KE^    = 


f2  RE3  y3  D3 


-i  0.333 


fl  K6 


3  r^3 

yl  D2h 


(16a) 


To  begin  the  design  procedure,  a  value  of  RE„  is  assumed. 
Use  the  f2  verses  RE  plot  for  side  2  to  determine  f2»   A  sample 
data  plot  from  reference  [2]  along  with  the  associated  geomet- 
ric parameters  describing  a  typical  side  2  flow  arrangement  is 
shown  in  Figure  A-II. 
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FIGURE   A- I I 


SAMPLE    DATA   PLOT   FOR   SURFACE    2 


[REF    2] 
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It  should  be  noted  that  the  friction  factor  verses  RE  plot 
for  flow  over  finned  tube  banks  can  be  closely  approximated  by 
a  straight  line  on  a  Log-Log  plot.   This  implies  that  f~  can  be 
accurately  represented  by  an  equation  of  the  form: 

Bf 
f2  =  Af [RE] 

Table  I  list  the  linearized  equations  of  f_  [and  j_]  for 
Figures  92-101  of  reference  [2] . 

This  simplification  allows  for  ease  in  programming  as  dem- 
onstrated in  the  suggested  design  procedure  flow  diagram  out- 
lined in  Table  II. 


For  exact  calculations,  f,  should  be  read  from  the  f,  ver- 
ses RE  plot  for  the  specific  tubular  surface  under  investiga- 
tion, if  such  data  is  available.   From  this  plot  for  surface  1, 
it  is  possible  to  determine  RE.,  and  f..  to  satisfy  equation  (16) 
or  (16a)  once  a  value  of  RE?  has  been  assumed. 

If,  after  assuming  a  value  of  RE,,,  the  required  value  of 
f,  and  RE,  lay  outside  the  range  of  data  for  surface  1,  a  new 
value  of  RE2  must  be  assumed. 

If  the  required  data  for  surface  1  is  not  available,  fric- 
tion factor  for  turbulent  flow  inside  smooth  tubes  is  well  cor- 
related and  given  by  Nikuradse  as: 


/i-  =4.0  Log   [RE  /f]  -  0. 4< 
/I  i  o 


(17) 
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The  transition  from  laminar  to  turbulent  flow  inside  smooth 
tubes  occurs  between  Reynolds  number  of  2300  and  4000. 

Excellent  linear  approximations  of  the  Nikuradse  relation 
are  given  by: 

-0.25 


f1    =  0.0791  RE 


2300  <  RE  <  20,000    (18) 


-0.20 


fl    -  0.0460  RE 


20,000  <  RE  <  10 


(19) 


In  a  crossflow  finned  tubular  heat  exchanger  the  dominant 
heat  transfer  resistance  will  probably  be  the  gas  outside  the 
tubes  flowing  over  the  fins  (side  2) .   It  is  for  this  reason 
that  high  accuracy  for  the  f  (and  j)  prediction  of  the  fluid  on 
side  1  is  not  critical. 

Note  that  equation  (16a) ,  using  the  linear  approximation 
for  f.  given  by  equation  (18)  may  be  rewritten  as: 

.3636 


RE.,  = 


f2  RE2  A   Dlh 


0.0791  K6  V\   D^h 


(20) 
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Having  satisfied  either  equation  (16)  using  f..  verses  RE.. 
data  for  the  specific  tubular  surface  under  investigation,  OR 
equation  (20)  using  the  specified  linear  approximation  for  tur- 
bulent flow  inside  smooth  tubes,  length  dimensions  X  and  Y  may 
be  calculated  from  the  relations  given  in  equation  (10) : 

3 

X  =    2   2h  I0    _2  (22) 

2  f2  RE*  v\ 

3 

AP.,  D.,  g   p, 

y  =    1       lh  1°      1  (23) 

2  f1    REJ  yj 


Length  Z  may  be  calculated  from  equations  (5)  and  (7) : 

4  W. 

Z  =   i (24) 

^l  RE1  yl  X 
The  number  of  tubes  can  be  determined  from: 


„  =  44-  =  ^  °lh  (25) 

REi  "i  \ 


As  was  the  case  for  friction  factor,  j,  should  be  read 
from  the  j,  verses  RE  plot  for  the  specific  tubular  surface  un- 
der investigation,  if  available. 
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In  the  event  that  the  actual  heat  transfer  data  for  side 
1  is  not  available,  correlations  can  be  used  without  introduc- 
ing significant  deviations  as  explained  earlier. 

The  Colburn  correlation  for  forced-convection,  turbulent 
flow  in  tubes  is  widely  accepted: 

-0.20 
j1  =  0.023  RE  (27) 

The  heat  transfer  coefficient  [h]  may  be  defined  as: 


jl  Col  yl  RE1 

h  =  _i — E± — ± L  (28) 

Pr'667  D 
rrl    ulh 


Or  using  the  Colburn  correlation  expressed  by  equation  (27) : 

0.80 
0.023  RE     C  .  u, 
h.  =  2±-±  (28a) 

-i667  °lh 


The  value  of  j.  should  be  read  from  the  j2  verses  RE  plot 
ice  2.   As  was  the  case  for  f.~, 
be  expressed  in  an  equation  of  the  form: 


for  surface  2.   As  was  the  case  for  f2,    the  Colburn  modulus  can 


j2  =  A.  RE   * 
Table  I  presents  the  linearized  equations  of  j2  for  Fig- 
ures 92-101  of  reference  [2] . 
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j2  Cd2  y2  ^2 
h   =  -t Ef £ £  (29) 

P  '667  D 
r     u2h 


Refering  to  Figure  A-III,  calculate  the  appropriate  value  of 
fin  and  surface  efficiency  for  side  2: 


- +  - (33) 


AU     Am   U,      Am   U2    A.J,   h,      Am   n^   h 


Tl   !       T2   2    ^1   X       T2   °2   2 


where : 


A    =  a   XYZ 

xl     x 


A^  =  a2  XYZ 


If  c.  <  c0,  C  .  =  Cw  c  .  /c    -  c./c, 
1     2    mm     l   mm   max    1  z 

T.    -  T 
in,     out, 

E  =  ± i  (34a) 

T.    -  T. 

inl     in2 


If  C2  <  Cl'  Cmin  "  C2'  Cmin/Cmax  =  C2/Cl 

T      -  T. 
out0    mn 

e  =  ± ±-  (34b) 

T.    -  T. 
inl     in2 
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FIGURE  A-III 


Compute  fin  efficiency  assuming  a 
Longitudinal  fin  of  Rectangular 
Profile: 

Tanh  m£ 


Vl 


~mT 


where: 


.50 


2  h. 


FT 
L  m 


(30) 


What  ocometric  shape  describes  thc 
fin  being  employed  on  the  tube 
surface? 


Radial  Fin  of  Rectangular 
Profile 


An 


correction 


CA   +   DA  ^i-J 


whera: 


C."  -.00259  -  .0084  In  — 
A  r 


D.«  -.02426  -  .4521  In  — 
A  r. 


[Ref   3   page   264] 


Papered  (Hyperbolic)  Rad- 
ial fin 


An  "0 

correction 


[Ref  4  pg. 3-116] 


■<-L 


n ^-corrected 
ncorrection 


Compute  Surface  n: 
—  N  -  n 


Rectangular  Profile  fin 

sheet 
[see   below  for  I   defn 


An 

correction 

Same  as  Radial  Fin  of 
Rectangular  Profile 
using  r  (staggered  bank) 


<  ■  M 


50 


Xd   Do 


I  =   r 


[Ref   4   pg. 3-116] 


(31) 


^-corrected 


(32) 
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also 

A  U 


NTU  =  (35) 

min 


If  side  2  is  in  the  form  of  rectangular  fin  sheets,  the 
gas  flowing  over  the  tubes  and  the  fluid  in  the  tubes  are  con- 
sidered UNMIXED.   Using  Figure  A- IV  with  C  .  /C    and  the  mag- 

3    3  min  max  ^ 

nitude  of  e  calculated  from  equation  (34a  or  34b)  read  the  mag- 
nitude of  NTU  required  to  produce  this  effectiveness.   Find 

[NTU.^,,.   , ]  [C  .  ]  and  compare  with  the  results  of  equation 
required   min 

(33).   If  the  equation  is  not  balanced,  a  different  value  of 
RE2  must  be  assumed  and  the  procedure  repeated. 

If  side  2  is  in  the  form  of  low  fins,  the  gas  flowing  over 
the  tubes  is  considered  MIXED  while  the  fluid  in  the  tubes  re- 
mains UNMIXED.  The  procedure  outlined  above  is  used  with  Figure 
A-V  to  find  the  NTU  required  to  produce  the  calculated  effect- 
iveness in  this  case.   An  alternative  to  using  Figure  A-V  is  to 
calculate  the  NTU  requirement  from  equation  (14)  or  (15)  of  ref- 
erence [2]  . 

Table  II  presents  the  entire  design  procedure  in  a  flow 
diagram. 
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FIGURE  A- IV 


[BOTH  FLUIDS  UNMIXED] 
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80 
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FIGURE  A-V 

[ONE  FLUID  MIXED;  ONE  FLUID  UNMIXED] 
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TABLE  I 

LINEARIZED  EQUATIONS  OF  THE  FORM  A  [RE]  B  FOR  FRICTION  FACTOR  & 
COLBURN  MODULUS  FOR  FINNED  TUBULAR  SURFACES  PRESENTED  IN  REFER- 

ENCE  [2] 


FRICTION  FACTOR  [f] 

COLBURN  MODULUS  [j] 

FIGURE 

Af 

Bf 

A. 

D 

B. 

92 

.1781 

-.1954 

.3311 

-.4804 

93 

.1862 

-.2016 

.2262 

-.4127 

94 

.2377 

-.2192 

.2254 

-.4027 

96 

.2646 

-.2559 

.1461 

-.3418 

97A 

.2319 

-.2454 

.1541 

-.3596 

97B 

.5416 

-.2676 

.2163 

-.3666 

98A 

.2263 

-.2707 

.0842 

-.3273 

98B 

.2905 

-.2531 

.0947 

-.3239 

98C 

.3850 

-.2462 

.1023 

-.3137 

98D 

.2173 

-.2590 

.0813 

-.3264 

98E 

.3020 

-.2607 

.0915 

-.3288 

99A 

.1542 

-.2478 

.1432 

-.3726 

99B 

.3599 

-.2563 

.2027 

-.3787 

100 

.1223 

-.2173 

.1650 

-.3986 

101 

.0966 

-.2270 

.0940 

-.3458 

NOTE;  These  equations  are  written  for  the  range: 

1000  <  RE  <  upper  limit  of  data  IAW  Ref  [2] 
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CROSSFLOW  FINNED  TUUUI.AR  HEAT  EXCHANGER  DESTGN  PROCEDURE 


OH 


Assume  a 
value  for  RE. 


INPUT 
Linearized  equation  for  f  and  j  of 

D 

the  form  AIRE)   for  the  enhanced 
surface:  (subscript  2) 

Af(e)'  Bf(e)'  Aj(e);  Bj(e) 

Geometric  Factors 


Dlh'  D2h' 


tv   e2, 


i*-n/ 


I.    t 

Design  Constraints: 
T 


V 


Ml'  Tini 


W2,  TinJ 


out1#  AP1'  Pinx 


out. 


2    in. 


Other i 


Properties  of  the  two  working 
fluids,  fin  material  [k]  and 
shape,  constants  used  to  cor- 
rect fin  efficiency 


AP, 


!1  W2 


Ari  °1  82  Wl 


Subscripts : 


1  Hotter  fluid  in 
Smooth  Tubes 

2  Colder  fluid 
flowinq  over  f in- 
ned  tube  banks 


RE, 


3+B 


Af(e)  ^2 


f(e)  i>3  D3 
v2   Dlh 


(.0791)  K6  u3  D3h 


0.364 


Assumption : 


f1=0.0791    RE~ 
J1=0.0230    RE~ 


0.25 


0.20 


4P2   D2h 


\>2 

2+B 


2  Af  (e)    ^2 


f  (e) 


TABLE    II 
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*t  -  VYZ 


+  .80 
0.027  REX     C    Vl 

1  jSTJ 

Pr    Dlh 


Calculate: 


hl  *T, 


NO 
Go  To  A  and  Assume 
another  value  for 
RE. 


4  W, 


ax  rex  ux  X 


flPl  Plh  «.  °2 
0.1582  ReJ-75  v\ 


Go  To   B  to  compute 
n 


AUreq'd-  NTU  Cmin 


DOES 


AU 


req'd 


hl  V 


h2  \  AT, 


A,  -  S2XYZ 


1+Bj(e) 
A.(e)  RE    3(e>Cp2  U2 

2 

£  .667 


Calculate : 


h2?n< 


YES 
End  of  design  Procedure 


Ill 


Compute  fin  efficiency  assuming  a 
Longitudinal  fin  of  Rectangular 
Profile: 

Tanh  md 


"i-L 

m£ 

where i 
in 

~2  h2n 

.51 

k  & 
_  m  . 

What  geometric  shape  describes  tht. 
fin  being  employed  on  the  tube 
surface? 


Radial  Pin  of  Rectangular 

Prof 

lie 

Ancorre 

ction 

cfi  -  DA  Inln^j.) 

whera i 

r 
CA«  -.00259  -  .0084  In  -^ 
A                   re 

D.«=  -.02426  -  .4521  In  — 
A                   re 

[Ref  3  page  264] 

Tapered  (Hyperbolic)  Rad- 
ial fin 


Rectangular  Profile  fin 

sheet 
[see   below  for  I   defn] 


Correction" 


(Ref  4  page  3-116] 


correction 

Same  as  Radial  Fin  of 
Rectangular  Profile 
using  r  (staggered  bank) 

■50   Xd  Do 
d   o 


l-M 


l   =  r  -  r 
e    o 

[Ref  4  page  3-116] 


nd-corrected 
'<-L     Ancorrection 


Compute  Surface  ns 
A,.. 


n  -  1 


K  I 


1  -  n 


^-corrected 
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TABLE  III 

SUGGESTED  PROGRAM  FOR  SIZING  CROSSFLOW  FINNED  TUBULAR  HEAT  EX- 
CHANGERS  FOR  A  GIVEN  JOB  I AW  TABLE  II* 


I.  Design  Constraints/Geometric  Constants  to  be  stored  in 
Memory  as  follows: 


Primary  Register 

Secondary  Register 

Other 

Quantity   Units 

Quantity   Units 

Quantity   Units 

0   AP2        LB/FT2 

0  CP2/P-667  BTU/LBM-F 

A  RE2 (assumed)   ND 

1   A((eJ      ND 

1  g  pm2    LBM/HR2-FT2 

B  0.00 

2   Bf(e)      ND 

2  W        LBM/HR 

C  0.00 

3  M2                   LB/HR-FT 

3  gx       FT2/FT3 

D  0.00 

4   Dlh        FT 

4  AP.      LB/FT2 

E  A^/AT          ND 

5   CA         ND 

A 

5  g  p      LBM/HR2 -FT2 

I  62          FT2/FT3 

6   DA         ND 

A 

c     /i    /p  •  6  6  7 

D  ^pl'  r    BTU/LBM-F 

7   Wx        LB/HR 

7  A . ,  %      ND 
D  (e) 

8   y         LB/HR-FT 

Cut2 
8   p2M2^     FT 

9   D2h       FT 

K  P-G676 

m  r 
9  Bj(e)      ND 

II.  OUTPUT 

RE,,  X,  Z,  Y,  An*.   w      i«^jj*  Tl/        *.  jr  [NTU  C  .  ]   ,  ,  .  , 
1'   '   '   '   'to  be  applied    ((-corrected'       min  cal'td 


*For  Hewlett-Packard  HP-67/97  Programmable  Calculator 
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HI  .       PROGRAM  FOR  SIZING  FINNED  TUBULAR  HEAT  EXCHANC-'RS  FOR  A  GIVEN 

JOB* 


[Table  III  continued] 

601 

»L8i.A 

21  11 

057 

X* 

53 

113 

FZS 

16-51 

169 

Si  CD 

35  14 

002 

?:i 

le'-51 

058 

i 

-24 

114 

RCL4 

36  64 

170 

PCS 

16-51 

003 

RCL1 

36  81 

059 

RCLn 

36  11 

115 

X 

-35 

171 

(?CL3 

36  83 

004 

RCL3 

3i  63 

060 

2 

82 

lib 

RCLS 

36  65 

172 

RCLC 

36  13 

005 

X 

-35 

0$  1 

RCL2 

36  £2 

117 

X 

-35 

173 

X 

-35 

886 

aa.4 

56  64 

862 

* 

-55 

119 

.  pa;/ 

-14 

174 

STOC 

55  13 

00  7 

t 

'24 

863 

V 

31 

•  119 

sac 

36  13 

175 

LS7X 

16-63 

608 

SCL5 

36   05  ' 

864 

T 

-24 

128 

X 

-35 

176 

RCLO 

36  14 

00? 

♦ 

-24 

865 

2 

02 

121 

S7CC 

35  13 

177 

X 

-35 

010 

RCL2 

36  82 

866 

+ 

-24 

122 

RCLrt 

36  11 

178 

RCLI 

36  46 

011 

♦ 

-24 

867 

RCL1 

3b  01 

123 

RCL9 

36  65 

179 

X 

-35 

012 

KL1 

36  46 

868 

T 

-24 

124 

1 

6! 

188 

STOO 

35  14 

813 

T 

-24 

869 

RCL9 

3b  a 

125 

* 

-55 

161 

CCLm 

36  11 

014 

f;s 

16-51 

070 

3 

03 

126 

V" 

31 

182 

RCL9 

36  65 

015 

»CL6 

36  08 

671 

v 

31 

127 

car 

36  87 

183 

1 

81 

816 

X 

-35 

872 

X 

-35 

128 

X 

-35 

184 

♦ 

-55 

017 

RCL7 

36  67 

873 

p:s 

16-51 

129 

2 

82 

185 

yx 

51 

018 

X 

-35 

674 

RCL1 

36  81 

136 

X 

-35 

186 

RCL" 

36  07 

019 

ST07 

35  er 

075 

X 

-35 

131 

RCLS 

36  es 

187 

X 

-35 

020 

•L0L3 

21  11 

876 

PRTX 

-14 

132 

X 

-35 

188 

RCL8 

36   60 

021 

RCL3 

36  63 

877 

STOC 

35  13 

133 

p:s 

16-51 

1*9 

X 

-35 

822 

RCL4 

36   04 

678 

4 

64 

134 

RCLS 

36   05 

198 

f:s 

16-51 

821 

X 

-35 

879 

RCL2 

36  61 

135 

i 

-24 

191 

RCL3 

36  85 

824 

RCLS 

3b  0ft 

680 

X 

-35 

136 

IX 

«•  54 

!<)/> 

X 

-35 

825 

"T 

-24 

601 

RCL3 

36  63 

137 

ST00 

35  14 

1<>3 

RCL9 

36  65 

826 

RCL9 

36  85 

882 

T 

-24 

138 

e" 

33 

194 

T 

-24 

827 

+ 

-24 

883 

(tCLB 

36  12 

139 

PCLO 

36  14 

195 

RCLO 

36   14 

828 

3 

63 

084 

T 

-c4 

140 

CHS 

-22 

196 

X 

-55  • 

829 

*•* 

31 

085 

prs 

16-51 

14J 

e* 

33 

197 

1/X 

c  * 

J- 

030 

RCLi 

36  81 

886 

RCLS 

36  65 

142 

- 

-45 

19$ 

sroo 

55  14 

831 

A 

-35 

067 

T 

-24 

143 

IS"/ 

16-63 

199 

RCLS 

36   12 

832 

RCLO 

36  11 

668 

RCLC 

36  13 

144 

SCO 

36   14 

260 

, 

-62 

033 

,  wrx 

-14 

899 

T 

•24 

145 

«* 

33 

201 

8 

86 

034 

3 

83 

898 

frtx 

-14 

146 

♦ 

-55 

202 

0 

66 

035 

RCL2 

36  62 

891 

rclc 

36  12 

147 

i 

-24 

263 

V 

31 

036 

* 

-55 

892 

X 

-35 

146 

BCLC 

36   14 

284 

g 

-62 

837 

Y* 

31 

893 

STOC 

35  13 

149 

T 

-24 

205 

0 

66 

838 

X 

-35 

894 

RCL4 

36  64 

158 

EN7f 

-*li 

206 

2 

62 

839 

. 

-62 

695 

3 

63 

151 

LH 

32 

207 

3 

85 

040 

0 

88 

696 

Y* 

31 

152 

RCL6 

36  66 

208 

X 

-55 

041 

7 

87 

897 

ftCLB 

36  12 

153 

CHS 

-22 

209 

RCL4 

36  64 

042 

9 

89 

898 

1 

81 

154 

X 

-35 

210 

T 

-24 

043 

1 

81 

899 

. 

-62 

155 

RCLS 

36  65 

211 

RCLS 

36  65 

044 

* 

-24 

188 

7 

67 

156 

♦ 

-55 

212 

X 

-35 

045 

RCL7 

36  67 

181 

5 

65 

157 

FR7X 

-14 

213 

prs 

16-51 

04$ 

^ 

-24 

102 

V 

31 

158 

- 

-45 

214 

RCL6 

36  66 

047 

2 

12 

183 

i 

-24 

159 

ilOSJ 

35  14 

215 

X 

-35 

848 

. 

-62 

164 

RCLS 

36  6i 

160 

PKTX 

-14 

216 

RCLC 

36   i .' 

849 

7 

87 

105 

X* 

53 

161 

1 

01 

217 

X 

-55 

050 

5 

65 

16b 

i 

-24 

162 

RCcO 

36  14 

218 

1/X 

e  - 

051 

1/X 

52 

107 

. 

-62 

163 

- 

-45 

219 

RCLD 

56   14 

052 

V 

3\ 

108 

1 

61 

164 

RCIE 

36  15 

228 

♦ 

-55 

053 

frtx 

-14 

109 

5 

05 

165 

X 

-35 

221 

1/X 

5f 

054 

sroB 

35  12 

110 

6 

62 

166 

1 

61 

222 

FRTX 

-14 

655 

rcl6 

36   00 

111 

2 

6i 

167 

XZ't 

-41 

223 

r'JS 

16-5i 

05b 

RCL5 

36   03 

112 

* 

--'- 

168 

- 

-45 

224 

(TN 

<"' 

*For  Hewlett-Packard  67/97  Programmable  Calculator 


IV.  OPERATION 
First  Run 


Select  a  Value  of 
A. 


and  store  in  register 
To  begin  the  program,  press  key  "A" 


Subsequent  Runs:   [As  necessary  until  the  calculated  and 

required  value  of  NTU  C  .   agree]  Select  a 

mm  J 

value  of  RE2  and  store  in  register  A.   To 
begin  the  program  press  key  "B" . 
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HEAT  EXCHANGER  DESIGN  PROBLEM 
APPENDIX  IV 
A.  Description 

Appendix  III  procedure  can  be  used  for  sizing  crossflow 
finned  tubular  heat  exchangers.   The  following  data,  taken  from 
reference  [2] ,  is  used  to  determine  the  required  heat-exchanger 
size  for  the  given  conditions. 

The  design  is  that  of  an  intercooler  for  a  5000  SHP  gas 
turbine  plant. 


DESIGN  CONSTRAINT 


Surface 


6 
6 

VAT 
I 

W 

T. 
in 

T 
out 

AP 

y 


SIDE  I 

SIDE  II 

Smooth  Tube 

Finned  Tube  11.32 

.737  SR  [Fig  106] 

.01224  FT 

.01152  FT 

N/A 

.00033  FT 

42.1  FT2/FT3 

270  FT2/FT3 

N/A 

.845 

N/A 

.225  IN 

400,000  LB/HR 

200,667  LB/HR 

60  °F 

260  °F 

82  °F 

80  °F 

1.554  PSI 

.306  PSI 

2.36  LB/HR-FT 

.0482  LB/HR-FT 

1.0  BTU/LB-°F 

.24  36  BTU/LB-°F 

62.38  LB/FT3 

.1755  LB/FT3 
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SIDE  I  SIDE  II 

PR  6.8  .70 

Km  N/A  100  BTU/HR-FT-°F 

F  .55  IN  N/A 


b  .79  IN  N/A 

.0560  IN2 
(From  Fig  106) 


Ac  .0560  IN2  Calculate 


B.  Calculations 


*1   =  Ac/Afr[4/V  *  [,056°  ^  ] 


[4] 


[.55  IN  X  .79  IN]  [.01224  FT] 
=  42.1  FT2/FT3  (1) 


32  =  270  FT2/FT3  (2) 


K1   =  [42.1]  [.01224/4]  =  .129  (3) 


K2  =  [270]  [.01152/4]  =  .778  (4) 


=  [400,000]  [.778]  m   12#038  (9) 

3    [200,667]  [.129] 


K   -  [1.554]  [.01224]  [62.38]  _  1017    Q  (11) 

4    [.3061]  [.01152]  [.1755] 
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K-  -  H2.038]2   =  >Q756  (13) 

[1917.8] 


K   =  [-0756]    =  4.333  X  10"5  (14) 

0    [12.038K 


CALCULATIONS  FOR  THE  FINAL  RE   SELECTED  WILL  BE  SHOWN; 


Select  RE2  =  5760 


From  the  data  of  Surface  2  [Figure  106]  read: 
f2  =  .021  j2  =  .0054 

In  order  to  be  strictly  consistant  with  the  design  presented  in 
Ref  [2]  ,  friction  factor  [f ..  ]  will  be  taken  from  Figure  29  of 
Ref  [2] ,  rather  than  employ  any  linear  approximations  suggested 
in  the  design  procedure: 

From  Figure  29  data  for  surface  1  and  Equation  (16) : 


RE   =  5073        f1   =  .0072 


=  R^y^  =  [5760]  [.0482]  =  24flQ0  LB/HR.FT2 
2     D2h         [.01152] 


m    [.306]  [.01152]  [144]  [32.2]  [3600]2  [.1755]      {22y 
[2]  [0.021]  [24,100]2 
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X  =  18.29  IN  (22) 


Y  =  5.34  FT  (23) 


Z  =  [400,0001/978,127]  =  2^Q    pT  (24) 

[.129]  [18.28/12] 


where  G  =  [5073]  [2-36^  =  978,127  LB/HR- 
[.01224] 


N  m    [18.29]  [24.96]  =  1Q51 
[0.55]  [[0.79] 


FT2 


(25) 


Once  again,  in  order  to  be  strictly  consistant  with  the  design 
presented  in  Ref  [2],  Colburn  Modulus  [j.]  will  be  taken  from 
Figure  28  of  Ref  [2] ,  rather  than  employ  any  linear  approxima- 
tions suggested  in  the  design  procedure: 

From  Figure  28  data  for  surface  1  @  RET  =  5073: 

jn  =  0.00502    (read  N  .  =  .0014  and  compute  j] 


h   =  [0.005]  [1.0]  [978,127]  [6.8]"*667 

=  1368  BTU/HR  FT2  °F  (28) 


h2  =  [0.0054]  [.2436]  [24,100]  [0.700]"*667  (29) 

=  40.21  BTU/HR  FT2  °F 
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m2  =  /(2)  (40.21)/  (.00033)  (100) 

=  49.3  FT"  (30) 


=  Tanh  [49.3]  [.225/12]  =  Q1Q8  (3Qa) 

6        [49.3]  [.225/12] 


Ancorrection  =  °  ^For  consistancy  with  Ref  (2)        (31) 

example] * 


nQ  =  1  -  .845[1  -  .788]  =  0.82  (32) 


+  L 


AUcalculated    t1368l  [711.55]     [40.21]  [0.82]  [4570.9] 


130,506  BTU/HR-°F 


(33) 


In  order  to  complete  the  design  procedure,  compute  AU    ired 
and  compare  with  equation  (33)  results: 


[AU]     .     -  NTU  C 
1   J  required        min 


Cmin  "  c«ii-/.44«  x-n  =  [200,667]  [.2436] 
mm    air  (side  II) 

=  48,883  BTU/HR-°F 


C     =  C      ,  .  ,   _.=  400,000  BTU/HR-°F 
max    water (side  I) 
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260  -  82    „  oft 
e  =  — )i Z±   =  o.89  (34b) 

260  -  60 


From  Figure  5  of  Reference  [2] ,  for  a  cross-flow  exchanger  with 

both  fluids  UNMIXED,  and  C  .  /C    =  0.122  with  e  =  0.89,  [or 
min  max  ' 

from  Table  4  of  Ref  [2]  with  proper  interpolation]  the  NTU    ,, 
is  found  to  be  2.67. 


NTUreq'dCmin  =  ^2  .  67]  [48  ,  883]  =  130,517  BTU/HR-°F     (35) 


Comparing  Equation  (33)  with  equation  (35) : 


C 
min 


NTU    ,  ,C  .   =  NTU   ,   ,  .  , 
req'd  mm      calculated 

130,517  BTU/HR-°F    =  130,506  BTU/HR-°F 


CLOSE  ENOUGH! 


Surface  2  is  a  rectangular  fin  sheet  and,  to  be  strictly 
correct,  a  correction  should  be  applied  to  permit  the  use 
of  equation  (30a)  as  outlined  in  the  design  procedure. [See 
Appendix  III,  Table  I]  In  order  that  the  results  of  this 
example  will  be  consistant  with  those  of  reference  [2] ,  a 
correction  of  0  will  be  assumed  pursuant  to  the  simplifying 
assumptions  as  stated  in  the  example  of  reference  [2]. 
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